Universal gap fluctuations in the superconductor proximity effect 
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Random-matrix theory is used to study the mesoscopic fluctuations of the excitation gap in a 
metal grain or quantum dot induced by the proximity to a superconductor. We propose that the 
probability distribution of the gap is a universal function in rescaled units. Our analytical prediction 
for the gap distribution agrees well with exact diagonalization of a model Hamiltonian. 
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A normal metal in the proximity of a superconduc- 
tor acquires characteristics that are typical of the su- 
perconducting state Jf|. One of those characteristics is 
that the quasiparticle density of states vanishes at the 
Fermi energy. This superconductor proximity effect is 
most pronounced in a confined geometry, such as a thin 
metal film or metal grain, or a semiconductor quantum 
dot. In that case, provided the scattering in the normal 
metal is chaotic, no excitations exist within an energy gap 
E g ~ h/r, where r is the typical time between collisions 
with the superconductor ^J-^| . 

If the coupling to the superconductor is weak (as for 
the point contact coupling of Fig. 1), the functional form 
of the density of states becomes independent of micro- 
scopic properties of the normal metal, such as the shape, 
dimensionality, or mean free path. Weak coupling means 
that t is much bigger than the time T erg needed for er- 
godic exploration of the phase space in the normal region. 
For a point contact with N 3> 1 propagating modes at 
the Fermi level e = 0, the density of states has a square 
root singularity at the excitation gap Ef, 
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For a ballistic point contact and in the absence of a mag- 
netic field, E g = cNS and A g = c'TV 1 / 3 ^, where c = 0.048 
and c' = 0.068 are numerical constants and S is the mean 
level spacing in the normal metal when it is decoupled 
from the superconductor. 

Equation (|l|) was obtained in a self-consistent diagram- 
matic perturbation theory that uses rS/h ~ N^ 1 as a 
small parameter. Such a mean-field theory provides a 
smoothed density of states for which energies can only be 
resolved on the scale o£h/r ~ NS, not on smaller energy 
scales, and is unable to deal with mesoscopic samplc- 
to-samplc fluctuations of the excitation gap. Mesoscopic 
fluctuations arise, e.g., upon varying the shape of a quan- 
tum dot or the impurity configuration in a metal grain. 
The lowest excited state e\ fluctuates from sample to 




FIG. 1. A quantum dot (N) connected to a superconduc- 
tor (S). The voltages on the gates V\ and V2 change the shape 
of the dot. Different values of the applied voltages create dif- 
ferent samples within the same ensemble. 

sample around the mean field value Eg, with a probabil- 
ity distribution P{s\)- It is the purpose of this paper to 
go beyond mean field theory and to study the mesoscopic 
fluctuations of the excitation spectrum close to Eg. Our 
main result is that the gap distribution P{e±) is a univer- 
sal function of the rescaled energy x = (£1 — E g )/A e , in 
a broad range \x\ <C N 2 / 3 . The Fermi level itself (e = 0) 
falls outside this range, which is why the universal gap 
distribution was not found in a recent related study [|9|. 
Our main findings are illustrated in Fig. 2. 

We first consider the gap distribution in the absence of 
a magnetic field, and then include a time-reversal sym- 
metry breaking magnetic field. Starting point of our cal- 
culation is the effective Hamiltonian [fbl 
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Here H is an M X M Hermitian matrix representing the 
Hamiltonian of the isolated quantum dot, and W is an 
M x N matrix that describes the coupling to the su- 
perconductor via an TV-mode point contact. For a bal- 
listic point contact, W mn = 7r _1 5 mn (M5) 1 / 2 @. The 
number M is sent to infinity at the end of the calcula- 
tion The effective Hamiltonian is a valid descrip- 
tion of the low-lying excitations if the Thouless energy 
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FIG. 2. Mean field and ensemble averaged density of 
states p m f and (p), together with the probability distribution 
P of the excitation gap, as a function of the rescaled energy 
x = (ei — Eg) / A s . These curves are the universal predictions 
of random-matrix theory. 



N8 is much smaller than the order parameter A of the 
bulk superconductor. In the absence of a magnetic held, 
the matrix H is symmetric. To describe an ensemble of 
chaotic quantum dots (or disordered metal grains), we 
take H from the Gaussian orthogonal ensemble (GOE) 
of random-matrix theory uM , 



V{H) oc exp 



A8 2 M 



TrH' 



(3) 



Calculation of the density of states of TL using perturba- 
tion theory in A -1 yields the result (Q) discussed in the 
introduction. Our problem is to go beyond perturbation 
theory and hnd the probability distribution P{e\) of the 
lowest positive eigenvalue ei of the Hamiltonian 

We have solved this problem numerically by exact di- 
agonalization of the effective Hamiltonian H. Before pre- 
senting these results, we first describe an entirely different 
approach, that leads to an analytical prediction for the 
gap distribution. We invoke the universality hypothesis 
of random-matrix theory, that the local spectral statis- 
tics of a chaotic system depends only on the symmetry 
properties of the Hamiltonian, and not on microscopic 
properties. This universality hypothesis has been proven 
for a broad class of Hamiltonians in the bulk of the spec- 
trum [Q, but is believed to be valid near the edge of 
the spectrum as well. A proof exists for so-called trace 
ensembles, having V(H) oc expj— tr/(JJ)], with / an ar- 
bitrary polynomial function |15| . 

The mean-field density of states near the edge can be 
written in the form 
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According to the universality hypothesis, the spectral 
statistics near the edge, in rescaled variables (e — b)/a 
depends only on the exponent p and on the symmetry 
index /?[/?= 1 (2) in the presence (absence) of time- 
reversal symmetry]. Generically, p is either 1/2 (soft 
edge) or —1/2 (hard edge). For our problem, we have 
P = 1, p = 1/2, a = 7T 2 / 3 A g , b = E g , cf. Eq. (@). The 
corresponding gap distribution is given by Jig] 
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Pi (a) = exp ( -\ I [q(x') + (x- x')q 2 (x')]dx' ) . (6) 
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The function q(x) is the solution of 

q "(x) = ~xq(x) + 2q 3 {x), 
Ai(- 



with asymptotic behavior q(x) 
[Ai(x) being the Airy function]. 

The distribution (H) is shown in Fig. 3 (solid curve). 
It is centered at a positive value of x = (e\ — E g )/A g , 
meaning that the average gapsize (e±) is about A g bigger 
than the mean-field gap E g . For small x there is a tail of 
the form 
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Non-universal corrections to the distribution (g) become 
important for energy differences |e — E g \ > E g , hence 
for |x| > A 2 / 3 . Since the width of the gap distribution 
is of order unity in the variable x, the probability to 
find a sample with an excitation gap in the non-universal 
regime is exponentially small. 

In order to verify our universality hypothesis, we com- 
pare Eq. (|^) with the results of an exact diagonalization 
of the Hamiltonian (Q). As one can see in Fig. 3, the nu- 
merical data are in good agreement with the analytical 
prediction. The small deviations can be attributed to the 
finiteness of N and M in the numerics. 

Let us now consider the effect of a weak magnetic field 
on the gap distribution. In the effective Hamiltonian, the 
presence of a magnetic field is described by replacing H 

by 



H(a) = H + iaA. 
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Here A is an M x M real antisymmetric matrix, whose 
off-diagonal elements have the same variance as those of 
H. The parameter a is proportional to the magnetic 
field, 
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where $ is the magnetic flux through the quantum dot, 
$o = h/e is the flux quantum, and r\ is a non- universal 
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FIG. 3. Probability distribution of the rescaled excitation 
gap x = (si — Eg)/A g . Data points follow from an exact di- 
agonalization of 10 4 realizations of the effective Hamiltonian 
(§) for different values of M and N (A: M = 400, N = 200, 
□: M = 600, N = 150, O M = 600, N = 80). The solid 
curve is the universal prediction (^) of random-matrix the- 
ory. The mean of the data points has been adjusted to fit the 
curve by applying an horizontal offset; no other fit parame- 
ters are involved. The inset shows the actual mean (a;) and 
root-mean-square value a of the data for M/N — 4 for differ- 
ent values of N, together with the random-matrix prediction 
for TV — > oo. These results are all in zero magnetic field. The 
dashed curve is the random-matrix theory prediction ( |l5| ) in 
the presence of a time-reversal-symmetry breaking magnetic 
field (/3 = 2). 



numerical constant |11|. The case a = corresponds 
to the GOE that we considered above; the case a = 1 
corresponds to the Gaussian unitary ensemble (GUE) of 
fully broken time-reversal symmetry. 

The effect of a magnetic field on the density of states 
in mean-field theory is known 0] . The square-root singu- 
larity ([j]) near the gap still holds, but the magnitude of 
the gap is reduced. The critical flux <!> c at which E g = 
and hence the proximity effect is fully suppressed is given 
by 
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This is a much larger flux than the flux $buik at which 
the spectral statistics in the bulk of the spectrum crosses 
over from GOE to GUE, which is given by |17[ 
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We will now argue that the characteristic flux $ c dgo 
for the spectral statistics at the edge of the spectrum 
is intermediate between <£> c and 'f'buik- We consider the 
effect of the magnetic field on the lowest eigenvalue E\ of 
H to second order in perturbation theory, 



8e\ 



^ Q . 2 l(iW>P 



A 



o 

-A 



(13) 



Since typically | 2) | 2 ~ MS 2 /tt 2 ande 2 -£i ~ A g , we 
see that the effect of level repulsion from the neighboring 
level £2 on the lowest level e\ becomes comparable to 
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The terms in Eq. (|13| ) with j> 1 give a uniform shift of 
all low-lying levels, and hence do not affect the fluctua- 
tions. For N 1 the flux scale (jlj) for breaking time- 
reversal symmetry at the edge of the spectrum is much 
smaller than the critical flux 4> c needed to suppress the 
proximity effect. What is needed is TV 2 / 3 <C N. This 
condition is difficult to satisfy in a numerical calculation. 
The analytical prediction for fully broken time-reversal 
symmetry is flfj| l 
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This curve is shown dashed in Fig. 3. The tail for small 
x is now given by 
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To make contact with Ref. || we briefly discuss the im- 
plications of our results for the ensemble averaged density 
of states (p(e)) in the sub-gap regime. The tail of P(x) 
for x < — 1 is the same as the tail of (p), cf. Fig. 2. We 
conclude that EM 
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over a broad range A g <C E g — e -C E g inside the mean- 
field gap. A different exponential decay (with a power 2 
instead of 3/2 in the exponent) was predicted recently by 
Beloborodov, Narozhny, and Aleiner Q, for the sub-gap 
density of states of an ensemble of superconducting grains 
in a weak magnetic field. Since the mean-field density 
of states in that problem is also of the form (Q), the 
universal GUE edge statistics should apply. The reason 
that the universal decay (|l^) was not obtained in Ref. ^ 
is that their theory applies to the non-universal energy 
range e -C E g near the Fermi level. To emphasize the 
significance of the universal energy range we note that 
the probability to have the lowest energy level in that 
range is larger than in the non-universal range by an 
exponentially large factor cx exp[(_E g /A g ) 3 / 2 ]. 

In conclusion, we have argued that the proximity effect 
in a mesoscopic system has a gap distribution which is 
universal once energy is measured in units of the energy 
scale A g oc (£ g <5 2 ) 1/3 defined fr om the mean-field density 
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Energy scale 


Flux scale 


Bulk statistics 
Edge statistics 
Gap size 
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TABLE I. Characteristic energy and magnetic flux scales 
for the spectral statistics in the bulk and at the edge of the 
spectrum and for the size of the gap. 



of states p(e) — \{e — _E g )/A|] 1 ' /2 /7r. This universal dis- 
tribution is the same as the distribution of the smallest 
eigenvalue of the Gaussian orthogonal or unitary ensem- 
bles from random-matrix theory, depending on whether 
time-reversal symmetry is broken or not. We have iden- 
tified the magnetic field scale for breaking time-reversal 
symmetry and verified our results by exact diagonaliza- 
tion of an effective Hamiltonian. Characteristic energy 
and magnetic field scales are summarized in Table 1. The 
universality of our prediction should offer ample oppor- 
tunities for experimental observation. 
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holds, but now with coefficients c and c' for the energy 
scales Eg and A g that weakly depend on M/N. For the 
comparison of the numerical data with the prediction of 
random-matrix theory we calculated E g and A g from the 
mean-field theory (Eq. (8) of Ref. [jjj) for the values of 
M/N used in the exact diagonalizations. No fit parame- 
ters are involved in this procedure. 
[13] M. L. Mehta, Random Matrices (Academic, New York, 
1991). 

[14] T. Guhr, A. Muller-Groeling, and H. A. Weidenmuller, 

Phys. Rep. 299, 189 (1998). 
[15] E. Kanzieper and V. Freilikher, Phys. Rev. Lett. 78, 3806 

(1997). 

[16] C. A. Tracy and H. Widom, Commun. Math. Phys. 159, 

151 (1994); 177, 727 (1996). 
[17] A. Pandey and M. L. Mehta, Commun. Math. Phys. 87, 

449 (1983). 

[18] The complete random matrix-theory prediction is 
{p(x)) = -xAi 2 (x) + [Ai'(x)] 2 + §6 p ,iM(x)[l - 
Ai(y)dy]. The = 1 result is plotted in Fig. 2. 



[1] M. Tinkham, Introduction to Superconductivity (Mc- 

Graw Hill, New York, 1995). 
[2] A. A. Golubov and M. Yu. Kupriyanov, Sov. Phys. JETP 

69, 805 (1989). 
[3] W. Belzig, C. Bruder, and G. Schon, Phys. Rev. B 54, 

9443 (1996). 

[4] J. A. Melsen, P. W. Brouwer, K. M. Frahm, and C. W. J. 

Beenakker, Europhys. Lett. 35, 7 (1997); Physica Scripta 

69, 223 (1997). 
[5] A. Lodder and Yu. V. Nazarov, Phys. Rev. B 58, 5783 

(1998). 

[6] F. Zhou, P. Charlat, B. Spivak, and B. Pannetier, J. Low 

Temp. Phys. 110, 841 (1998). 
[7] W. fhra, M. Lead beater, J. L. Vega, and K. Richter, 

|cond-mat/9909100[ 
[8] In a quantum dot or metal grain of size R, with Fermi 

velocity vf and mean free path £, one has l/r CIg ~ 

v F R- 2 mm(£,R). 
[9] I. S. Beloborodov, B. N. Narozhny, and I. L. Aleiner, 

|cond-mat/99122l4 

[10] K. M. Frahm, P. W. Brouwer, J. A. Melsen, and C. W. 
J. Beenakker, Phys. Rev. Lett. 76, 2981 (1996). 

[11] C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

[12] In our numerical computations it was necessary to choose 
the ratio M/N relatively small in order to achieve the 
limit JV > 1 needed for universality of the gap distri- 
bution. For finite M/N , the mean-field result (Q) still 



4 



